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Abstract 

Some topics from recent progresses in lattice QCD are reviewed. 
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1. Lattice QCD at T — 

1.1. Full QCD simulations near or at the physical quark mass 

Due to large computational costs, u, d quark mass in current lattice QCD simulations 
is heavier than the physical value. One then has to extrapolate results obtained at heavier 
quark mass to those at the physical point, using some formula such as chiral perturbation 
theory (ChPT). It is not clear, however, whether ChPT can be applied at rather heavy 
quark masses employed so far in lattice QCD simulations. This situation is gradually 
changing, thanks to improvements for both computers and algorithms. Indeed the first 
important message from the annual lattice conference of this year, lattice 2008, is that 
full QCD simulations near the physical light quark mass becomes possible. 

In tabl dl.ll parameters such as number of flavors, the lattice spacing a , the spatial 
extension of the lattice L, the minimum pion mass m™ 111, and L x m™ ln< , are listed for 
recent large scale full QCD simulations. In the left half of the list, the conventional quark 
action such as the Wilson quark action or the staggered quark action is employed for 
simulations, while the chirally symmetric quark action such as the overlap quark action 
or the domain-wall quark action is used in the right half. Let me point out here that 
the minimum pion mass in the full QCD simulations now becomes below 200 MeV for 
conventional actions and close to 300 MeV for chirally symmetric actions, though finite 
size effects could be sizable in some cases where Lm^ < 3. 
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Table 1 

(Incomplete) list of recent full QCD simulations. 
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Results for hadron spectrum from one of the most extensive calculations have been 
reported in Ref.[3J, where both chiral and continuum extrapolations have been performed. 
In this calculation, the minimum pion mass m™ m ' = 190 MeV is light enough, and m^L > 
4 is always satisfied. In Fig[Tl the hadron spectrum obtained in the continuum limit is 
compared with experiment. Here the overall scale, the light quark mass and the strange 
quark mass are fixed by m=, and itik- As can be seen from the figure, the agreement 
between lattice QCD and experiment is excellent. 

It is also reported recently that the pion mass can be smaller than 160 MeV in the 2+1 
flavor full QCD simulation at a — 0.09 fmpQ. In Fig. [21 the pion mass squired divided 
by the light quark mass, m^/m u( i, is plotted as a function of the light quark mass, 
where black solid circles are previous results from CP-PACS/JLQCD collaborations [T2"] 
obtained at heavier quark masses while red solid circles are newt ones from PACS- 
CS collaboration[T|. As quark mass decreases, this ratio increasingly deviates from the 
straight line obtained from the fit with black circles, suggesting an existence of the chiral 
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Fig. 1. Hadron spectrum in the continuum limit, , ran and m= are used as inputs to fix free parameters 
of QCD. 
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Fig. 2. m% /m u d as a function of m u d at a = 0.09 fm and L = 2.9 fm. Red circles correspond to m-x = 702, 
570, 412, 296, 156 MeV from right to left. The vertical dashed line is the physical point. 

logarithm in the ratio. Note however that Lm™ m - = 2.3 in this calculation, so that the 
finite size effect might be sizable at the lightest quark mass. To reduce the magnitude 
of the finite size effect less than a % level, gauge configurations at L — 5.8 fm with 
m v ~ 140 MeV corresponding m^L ~ 4.1 are currently accumulated by this group. This 
simulation indeed will become the real QCD calculation. 

1.2. Lattice QCD and chiral perturbation theories 

The second message of this year from lattice QCD is that NLO SU(2) chiral perturba- 
tion theory(ChPT) seems to work at < 450 MeV while NLO SU(3) ChPT seems to 
fail at the physical strange quark mass. The first positive statement can be established by 
both 2-flavor full QCD simulations with the small enough quark mass and use of chirally 
symmetric quark such as overlap or domain-wall quark. The second negative statement 
has been concluded from analyses in 2+1 flavor full QCD simulations. As a consequence 
SU(2) ChPT analysis is applied even to data in 2+1 flavor QCD. 

In Figj3]mJ/m g (top-left) and f„ (bottom- left) are plotted as function of the quark mass 
m 2 with the 2-flavor overlap fermion at a = 0.12 fm[TT]. The lightest point corresponds 
to rriyr = 290 MeV while the heaviest to 720 MeV. The lightest 3 points at < 450 
MeV are fitted by NLO SU(2) ChPT formula given by 



where M = m q / '(AttJ) 2 (green dash-dotted line), M = to^/(47t/) 2 (blue dashed line) or 
M = m 2 /(47r/ 7r ) 2 (red solid line). Differences among these three choices are higher order 
of ChPT in eq.([T]). The fact that all three fits work reasonably well at m w < 450 MeV 
establishes the validity of the NLO SU(2) ChPT at m n < 450 MeV. It is also shown 
in Ref. [13] that all data points up to the one at = 720 MeV can be fitted by the 
NNLO formula of the ChPT with M = m 2 /(47r/) 2 , whose NLO fit well describes data 
even beyond the fitted range. The convergence behaviour of the chiral expansion is also 
plotted for m^/m q (top-right) and /„• (bottom-right) in Fig[3j which shows, for example, 
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Fig. 3. Left: /m q (top) and (bottom) as a function of m 2 at a = 0.12 fm. Right: NNLO chiral fit 
with M = m^/(47r/) 2 using all data points for m 2 /2Bm q — l(top) and fn/ f — l(bottom). Solid curves 
are the full NNLO fit, while dashed ones represent the NLO contributions. 

that the NLO contribution is a -10%( +28%) level for ml/m q at ~ 500 MeV, 
while the NNLO correction gives about +3%(18%). Although the chiral expansion is at 
least convergent for both quantities, the NNLO contribution for /„. becomes significant 
already at ~ 500 MeV, corresponding to the K meson mass region for the 3 flavor 
case. 

In the 2+1 flavor full QCD simulations with the domain- wall quark, it is reported[8] 
that the NLO SU(3) (partially quenched) ChPT can not explain the quark mass depen- 
dence of mass and decay constant in the pseudo-scalar meson sector around the physical 
strange quark mass region. Instead the formula in the NLO SU(2) ChPT plus the heavy 
strange quark theory works better even for these data in the 2+1 flavor full QCD. This 
fact is consistent with the conclusion in the previous paragraph that the NNLO correction 
in the SU(2) ChPT can not be neglected at ~ 500 MeV for these quantities. 

1.3. Topology 

The 3rd message of this year from lattice QCD is that it becomes difficult to change 
topological charge in full QCD simulations at the lighter quark mass and/or near the 
continuum limit . One of possible solutions to this problem recently investigated is QCD 
with fixed topological charge. |14ll5j . For example, the two point function between fla- 
vor singlet pseudo-scalar densities P(x) at fixed topological charge Q behaves at large 
separation as 

^J m P( X )mP(0)) Q = + ^ (tf - J*.) + (i,) , (2) 

where V is the space-time volume, m is the quark mass, (0)q is the expectation values 
at fixed Q, Xt — (Q 2 )/V an d C4 = (Q A ) C /V are the topological susceptibility and the 
4-th cumulant defined in the ordinary QCD at — 0. This formula tells us that one can 
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Fig. 4. Topological susceptibility xt as a function of sea quark mass m q for Nf = 2 (circles) and 
Nf =2 + 1 (triangles) at Q = 0. The solid lines are fits with the form of ChPT. 

extract informations about QCD at 9 = such as Xt or c 4 from behaviours of correlation 
functions calculated at fixed Q. This strategy is suitable for lattice QCD since it is 
easier to fix the topological charge Q than to change it correctly during simulations, in 
particular, at small quark mass and/or near the continuum limit. 

The JLQCD has taken this strategy for the dynamical overlap project, where the fixed 
topology not only guarantees the locality of the overlap Dirac operator but also makes full 
QCD simulations with overlap quarks much faster. In FigHJ the topological susceptibility 
Xt extracted from eq.([2]) at Q — is plotted as a function of the light quark mass m q 
for 2 flavor [TB] and 2+1 flavor [T7j QCD. The results show a clear linear dependence on 
the sea quark mass expected from the ChPT as Xt — m q 'E/2 for 2 flavor and Xt = 
m q m s Yi/ (2m s + m q ) for 2+1 flavor, where £ is the chiral condensate at massless limit 
and m s is the (fixed) strange quark mass. The fits to data using this formula, denoted 
solid lines in the figure, yield, after a renormalization, E(2 GeV) = [245(5) (10) MeV] 3 
(N f = 2) and E(2 GeV) = [240(5) (2) MeV] 3 (N f = 3). 

2. Lattice QCD at finite temperature and density 

2.1. Isentropic Equation of states at finite density 

It has been shown that the zero- viscosity (perfect) hydro calculations explain the RHIC 
experimental results well. This indicates that in heavy-ion collision, after thermalization, 
the system expands and cools with constant entropy. In this case, the equation of state 
(EoS) along the line of the constant entropy per baryon number, called the isentropic 
EoS, becomes more relevant than the ordinary EoS as a function of temperature at the 
fixed chemical potential /i. The entropy per baryon number, S/Nb is roughly equal to 
300 at RHIC, 45 at SPS and 30 at AGS. The lattice QCD at non-zero chemical potential, 
however, is notoriously difficult to be simulated. So far there is no established method 
to simulate the full QCD at arbitrary value of the chemical potential /i. If the chemical 
potential is small, the Taylor expansion in terms of /i works well 1 18]. 

The MILC Collaboration has calculated the isentropic EoS using the Taylor expansion 
method in 2+1 flavor QCD with the staggered quarks at ~ 220 MeV |19l20] , In 
FigJSjLeft) their result for p/T 4 where p denotes the pressure is plotted as a function 
of T along lines of the constant S/Nb — 30, 45, 300 and oo, at JV( = 4 (open symbols) 
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Fig. 5. Left: Pressure p/T 4 as a function of temperature T along the lines of constant entropy per baryon 
number in 2+1 flavor QCD [20l . Right: The ratio of pressure and energy density p/e as a function of e 
along the lines of constant entropy per baryon number 2+1 flavor QCD |2T1 . 



and N t = 6 (solid symbols) where N t is the number of the lattice sites in the temporal 
direction which controls the lattice spacing as a = l/(N t T). The result shows that the 
lattice discretization errors are small for this quantity and p/T becomes larger at fixed 
T as S/Nb decreases. 

In Figj5^Right) the result for p/e, obtained by RBC-Bielefeld Collaborations in 2+1 
flavor QCD with the staggered quarks at ~ 220 MeV using the Taylor expansion 
method [21]. is plotted as a function of e on lines of the constant S/Nb = 30, 45, 300 at 
Nt = 4 (filled) and N t = 6 (open), where e is the energy density. Form this quantities 
the velocity of sound c s , which is important for the hydrodynamics calculation in heavy 
ion collisions, can be extracted as c 2 s = 4^ = e d ^/ e ' + f . Although discretization errors 
are observed, the dependence on S/Nb seems small in this ratio. 



2.2. QCD critical point at finite density 

The finite temperature phase "transition" in 2+1 flavor QCD at the physical light 
(u,d) quark mass and the physical strange quark mass seems a cross-over at zero or small 
chemical potential, while the phase transition is expected to become of first order at large 
chemical potential. Between the cross-over region and the first order region, there must 
exist an unique second order phase transition point, which is called the critical end point 
of QCD. Important questions to be answered are as follows. Does the critical end point 
indeed exist ? If so at which temperature and density does it occur ? To answer these 
questions, lattice QCD simulations at large /i are required. Unfortunately the Taylor 
expansion method mentioned before does not work at such large /x. The reweighting 
method, which looked promising on small volumes, may not be reliable to answer these 
questions 22 . A new method appears [23]. but we have to wait for a while, in order to see 
whether the method indeed works for this problem. One established method to investigate 
the existence of the critical end point is the imaginary chemical potential method, which 
I will discuss in this subsection. 

At a given chemical potential \i in 2+1 flavor QCD with the fixed strange quark mass 
m s , the second order phase transition appears at m q = m c (/i) where m q is the light 
(u,d) quark mass. Furthermore the phase transition is of first order at m q < m c {p) 
while it becomes a cross-over at m q > m c {jj). The statement that the phase transition 
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Fig. 6. Left: The central potential for the 2N state with anti-periodic boundary condition at E ~ 50 
MeV(red bars) and periodic boundary condition (blue crosses) at E ~ in quenched QCD. Right: The 
central potentials for ^Sq (red) and 3 Si (blue) 1 = 1 NS states in quenched QCD. 



in 2+1 flavor QCD at /i = is a cross-over at physical m q and m s correspond to the 
fact that mP hys > m c (0) at m s = mP hys where m^ hy and m^ hys are physical light and 
strange quark masses, respectively. If m c (p) becomes larger than m c (0) at larger /x, it is 
expected that m c (/i) at m^ hys becomes equal to m^ hys at some \x = fi c ; m^ hys = m c (/z c ). 
This fi c corresponds to the value of the chemical potential at the critical end point of 
QCD. If this scenario is correct, it is expected that m c (p) increases from m c (0), so that 

d2 ™^ > o at (i = 0. (Note that = at ^ = 0.) 

dpi 2 - o/i 

The imaginary chemical potential method avoids the complex phase problem of lattice 
QCD with non-zero /i by taking /i = ifij, and lattice QCD with real /ir can be easily 
simulated. If m c (n) can be analytically continued to fj, = ifij, we have 

! ^-i+£*(&f->+E*<-i>'(&f- < 3 » 

v ' k=l k=l 

Therefore Ck can be extracted from m c (ni). The latest result [2^] in 3 flavor QCD with 
vrtq — m s at Nt = 4 yields 

which strongly indicates that the expected critical endpoint is absent at small fi/T. Of 
course it is necessary to check this unexpected behaviour by other methods for the definite 
conclusion. 



3. Potential between baryons 

The force between nucleons (nuclear force) play a central role in nuclear physics. In 
a recent paper [25] , which has received general recoginition 26], the potential between 
nucleons has been calculated from lattice QCD in the quenched approximation. The 
result qualitatively reproduces all the features of phenomenological potentials that the 
force at medium to long range is attractive while it becomes repulsive at short distance 
forming a characteristic repulsive core. 
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The potential between nucleons, V(x), is extracted from the wave function y>£;(x) 
through the Schrodinger equation as 



[E - %g(x) 



V 2 



¥> B (x) = (0\N(x,0)N(0,0)\2N,E) 



(5) 



Ho 



'2/j 



where N(x.,i) is a nucleon interpolating field, \2N,E) is a 2 nucleon (2N) eigenstate in 
QCD with energy E and fi is a reduced mass of 2N system (/i = itin/2). In Rcf.|25j, _E ~ 
has been used to extract the potential, but the potential from this definition may depend 
on the energy E of the 2N state. In Fig[6] (Left), the potentials obtained at E ~ O(blue) 
and at E ~ 50 MeV(red) are compared in quenched QCD [28] . The potentials are almost 
identical between two energies, tough statistical fluctuations becomes sizable at large r 
in the case of E ~ 50 MeV, partly due to the contamination from higher energy excited 
states. It seems that the energy dependence of the potential is weak at this energy range. 

An important applications of this method is to calculate the potential between nu- 
cleon and hyperon or between two hyperons, which can not be extracted directly from 
scattering experiments. In Fig[6] (Right), the potentials between N and S in the 1=1 
channel for 1 So (red) and 3 Si (blue) states in quenched QCD are presented [27] • While 
overall features are similar to those of nuclear potentials, stronger spin dependence of 
the potentials is observed in this case. 
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